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ABSTRACT 

We present new VLBI images of supernova 1986J, taken at 5, 8.4 and 22 GHz between t = 22 
to 25 yr after the explosion. The shell expands oc £°- 69±0 03 . We estimate the progenitor's mass- 
loss rate at (4 ~ 10) x 10~ 5 M® yr -1 (for v w — 10 km s _1 ). Two bright spots are seen in the 
images. The first, in the northeast, is now fading. The second, very near the center of the projected 
shell and unique to SN 1986J, is still brightening relative to the shell, and now dominates the VLBI 
images. It is marginally resolved at 22 GHz (diameter ^0.3 mas; ~ 5xl0 16 cm at 10 Mpc). The 
integrated VLA spectrum of SN 1986J shows an inversion point and a high-frequency turnover, both 
progressing downward in frequency and due to the central bright spot. The optically-thin spectral 
index of the central bright spot is indistinguishable from that of the shell. The small proper motion 
of 1500 ± 1500 km s _1 of the central bright spot is consistent with our previous interpretation of it as 
being associated with the expected black- hole or neutron-star remnant. Now, an alternate scenario 
seems also plausible, where the central bright spot, like the northeast one, results when the shock front 
impacts on a condensation within the circumstellar medium (CSM). The condensation would have 
to be so dense as to be opaque at cm wavelengths (~ 10 3 x denser than the average corresponding 
CSM) and fortuitously close to the center of the projected shell. We include a movie of the evolution 
of SN 1986J at 5 GHz from t = to 25 yr. 

Subject headings: supernovae: individual (SN 1986J) — radio continuum: supernovae 



1. INTRODUCTION 

Supernova 1986J was one of the most radio luminous 
supernovae ever observed. Its long-lasting and strong 
radio emission and its relative nearness makes it one of 
the few supernovae for which it is possible to produce 
detailed images with very-long-baseline interferometry 
(VLBI), in addition to being one of the few supernovae 
still detectable more than t = 20 years after the explo- 
sion. 

SN 1986 J was first discov ered in the radio , some 
time after the explosion ()van Gorkom et all 119861 : 
IRupen et al.l ll987). The best estimate of the explosion 
epoch i s 1983.2 ± 1.1 (t = 0, Bietenholz et al. 2002 
see a lso IRupen et aHll987t IChevalier|[l987t IWeiler et ail 
11990ft . It occurred in the nearby galaxy NGC 891 
whose distance is 
Ferrarese et al.l 
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We will adopt the round value of 
10 Mpc throughout this paper. 

Early optical spectra showed that the supernova was 
unusual, but based on its promin ent Ha lines it wa s 
classifi ed as a Type Iln s upern ova ()Rupen et al.l fl987). 
Later, iLeibundgut et al.l (|1991ft suggested that it might 
be of Typ e lb. VLBI obser vations started soon after 
discovery (Bartcl e t al.l Il987h . and an image, the first 
of any optically identifie d supernova, was obtained soon 
after (|Bartel et al.l 119911) . and showed a heavily modu- 
lated shell structure with protrusions. VLBI observa- 
tions at subsequent epochs up to 1999 led to a series of 
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images showing the expansion of the shell, and allowing 
the expansion velocity and deceleration to be measured 
(jBietenholz et alj|2002|) . 

The evolution of the radio spectrum of a supernova is 
in most cases predictable. In particular, once supernovae 
have become optically thin after a few months or years, 
they usually display power-law radio spectra with spec- 
tral indices, a, in the range of —0.8 to —0.5 (S oc v a ). 
SN 1986J also showed such a spectrum before 1998, but 
a remarkable change occurred after that. An inversion 
appeared in the spectrum, with the brightness increasing 
with increasing frequency above ~ 10 GHz, and a high - 
frequency turnover at ^20 GHz ( Bietenh olz et al.| [2002') . 

This spectral inversion was associated with a bright 
spot in the (pro jected) center of the ex panding shell, 
as we showed in IBietenholz et all (|2004aD using phase- 
referenced multi-frequency VLBI imaging. At that time 
(late 2002) the bright spot was clearly present in the 
15 GHz image, but not discernible in the 5 GHz one. 
We report here on new VLBI and VLA observations of 
SN 1986J which show the further evolution of this unique 
object. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. VLBI Measurements 

The VLBI observations of SN 1986 J were made be- 
tween 2005 and 2008 with global VLBI arrays of 13 to 
18 antennas and with durations of 12 to 15 hours. The 
observations of SN 1986J were interleaved with ones of 
3C 66A, only 40' away on the sky, which we used as a 
phase-reference source. Details of the four observing runs 
are given in Table[Q The declination of +42° of SN 1986J 
enabled us to obtain dense and only moderately elliptical 
u-v coverage. As usual, a hydrogen maser was used as a 
time and frequency standard at each telescope. The data 
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Table 1 

VLBI Observations of SN 1986J 



Date 


Frequency 


Antennas a 


Total time 


Recording rate 




(GHz) 




(hr) 


(Mbits s" 1 ) 


2005 Apr 25 


22 


VLBA, Ef, Gb, Y27 


12 


256 


2005 Oct 24 


5 


VLBA, Ef, Gb, Y27, Jb, On, Wb, Tr 


12 


256 


2006 Dec 3 


8 


VLBA, Ef, Gb, Y27 


15 


512 


2006 Dec 10 


22 


VLBA, Ef, Gb, Y27 


15 


512 


2008 Oct 26 


5 


VLBA, Ef, GB, Y27, Jb, Mc, Nt, Tr, Wb 


18 


512 



a VLBA, ton 25 m dishes of the NRAO Very Long Baseline Array; Ef, 100 m, MPIfR, Effelsberg, 
Germany; Gb, 105 m, NRAO, Green Bank, WV, USA; Y27, equivalent diameter 130 m, NRAO, near So- 
corro, NM, USA; Jb, 76 m, Jodroll Bank, UK; Mc, 32 m, IdR-CNR, Medicina, Italy; Nt, 32 m, IdR-CNR, 
Noto, Italy; On, 20 m, Onsala Space Observatory, Sweden; Tr, 32 m, Torun, Poland; Wb, equivalent diam- 
eter 94 m, Wcstcrbork, the Netherlands. 



were recorded with either the VLBA or the MKIV VLBI 
system with sampling rates of 256 or 512 Mbits per sec- 
ond. We observed at frequencies of 5.0, 8.4 and 22 GHz. 
The data were correlated with the NRAC0 VLBA pro- 
cessor in Socorro, NM (USA), and the analysis was car- 
ried out with NRAO's Astronomical Image Processing 
System (AIPS). The initial flux density calibration was 
done through measurements of the system temperature 
at each telescope, and improved through self-calibration 
of the 3C 66A data. 

We phase-referenced to 3C 66A for all of the present 
observing runs. Our positions in this paper are given rel- 
ative to an assumed position of RA = 02 h 22 m 39!611500, 
decl. = 43° 02' 07'.'79884 taken fr om the Interna tional 
Celestial Reference Frame, ICRF (jFev et al.ll200l . for 
the brightness peak of 3C 66A. 

2.2. VLA Observations 

In addition to the VLBI observation, we also obtained 
VLA observations to measure SN 1986J's total flux den- 
sity at a range of different frequencies, both during our 
VLBI observations and at other times. We also reduced 
a number of data sets from the VLA archives in order 
to obtain additional flux density measurements. The 
epochs and the frequencies observed are listed in Table [21 
These observations allowed us to monitor the evolution 
of SN 1986J's integrated radio spectrum. 

The VLA data were reduced following standard pro- 
cedures, with the flux density scale calibrated by us- 
ing observations of the standard flu x density cal i brator s 
(3C 286 and 3C 48) on the scale of iBaars et all (fl977i ). 
An atmospheric opacity correction using mean zenith 
opacities was applied, and NGC 891/SN 1986J was self- 
calibrated in phase to the extent permitted by the signal- 
to-noise ratio for each epoch and frequency. At lower 
resolutions the supernova cannot be reliably separated 
from the diffuse emission from the galaxy. Therefore, we 
cite only flux densities obtained at FWHM resolution^ 
of < 7". For the 0.33 GHz measurements we incorpo- 
rate an additional uncertainty to reflect the difficulty of 
accurate galaxy subtraction. 

We note also that the gradual introduction of EVLA 

4 The National Radio Astronomy Observatory, NRAO, is a facil- 
ity of the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 

5 We take the geometric mean of the major and minor axes 
of the elliptical Gaussian fit to the dirty beam as our measure of 
resolution. 



antennas into the VLA array might cause some prob- 
lems in the flux-density calibration. Since the band- 
pass of the EVLA antennas differs from that of the 
older VLA-antennas, non-closing offsets are introduced 
for continuum observations such as ours, with the gain 
of an EVLA-EVLA baseline being typically several per- 
cent higher than that of a VLA- VLA baseline. Such 
non-closing offsets cannot be calibrated out using only 
standard procedures, which assume that the gains are 
dependent only on the antennas, and not on the baseline. 
These closure errors might compromise the accuracy of 
the flux density calibration. To determine whether this 
effect is significant, we carried out, for several epochs 
and frequencies, a more elaborate flux-density calibra- 
tion scheme using baseline-dependent factors calculated 
by the AIPS task BLCAL. Although the more elaborate 
calibration led to an improvement in the dynamic range, 
in no case was the normal flux-density scale in error by 
more than 1%. Since we include in all cases a flux-density 
calibration uncertainty of > 5%, we can say that any er- 
rors introduced by the use of the EVLA antennas are 
well within our stated uncertainties. 

In Table [5] we list the total flux density measurements 
that we obtained from VLA observations. Note that 
som e of these values have been previously published 
(see IBietenholz fc Bartel l2008al IBietenholz et al.l 120051 
I2004air but we repeat them here for the convenience of 
the reader and to allow a better interpretation of the 
evolution of SN 1986J's radio emission. 

We plot a selection of these flux densities in two ways: 
first, in Figure[T] we plot the evolution of the flux-density 
as a function of time at several different frequencies. Sec- 
ond, in Figure [2] we plot the radio spectrum as a func- 
tion of time at a number of epochs. We focus in this 
paper on the flux-density evolution in the last decade, 
for lightcurves covering the ea r lier p art o f SN 1986J's 
evolut ion, see IBietenholz et al.l (|2002f) and iWeiler et al.l 

FigurcQ]shows that the flux density at all our observing 
frequencies is generally decreasing. In terms of a power- 
law decay with S u oc t" , a weighted least-squares fit gives 
an average value of f3 — —2.2 over the last decade and 
over the frequencies of 1.5, 5.0, 8.4, 15 and 22 GHz, with 
the fastest decay seen at 1.5 GHz with (3 = —3.3 and the 
slowest at 15 GHz with f3 = —1.0. 



SN 1986J VLBI 



3 



Table 2 

Flux Densities of SN 1986J from VLA Observations 





Date 


n on f ' T T 
0.32 GHz 


1.43 GHz 


Flux densities' 1 
4.9 GHz b 8.4 GHz c 


14.9 GHz 


OO C** TU d 


43.3 GHz 


1995 


Jun 22 c 




25.6 ± 1.3 


14.2 ±0.8 




9.1 ±0.7 


8.0 ±0.8 






1996 


Oct 25 e 






11.0 ±0.6 


8.7 ±0.5 










1997 


Jan 31 e 








8.2 ±0.5 










1998 


Feb 9 e 










5.65 ± 0.69 


9.3 ±1.5 






1998 


Jun 5 




13.7 ± 0.9 f 


8.5 ± 0.5 


7.2 ± 0.4 










1999 


Feb 2 






7.3 ± 1.5 


6.1 ± 0.4 


6.1 ± 0.1 


9.5 ± 1.0 






1999 


Jun 13 e 




11.69 ± 0.62 




5.68 ± 0.30 


4.92 ± 0.52 








2000 


Oct 19° 




8.97 ± 0.47 


5.43 ± 0.29 


4.76 ± 0.26 


5.88 ± 0.42 


5.67 ± 0.65 






2001 


Jan 25° 




9.02 ± 0.22 




4.75 ± 0.36 










2002 


Jan 27° 




7.68 ± 0.46 


4.15 ± 0.26 


4.13 ± 0.28 


5.02 ± 0.54 


4.22 ± 0.44 






zuuz 


Moir 9^ 




7 1/1 4- n 70 
1 . 1U zn u. / u 


a 9n 4- n qo 
4.ZU ± U.oU 


q en 4- n on 
o.ou ± U.oU 


A 70 4- ^0 


o.uu zn u.du 




q+2.1 
o,3 -0.9 


2002 


Nov 11 






3.86 ± 0.28 


A A c 1 n DC 

4.1b ± 0.25 




£ o'j _l_ n c;n 
0.0/ ± U.OU 




2003 


Jan 8 






-3.75 ± U.45 


o.7o ± U.oU 


1 1 n A A 

4.25 ± 0.44 


A A 1 1 n A 

4.41 ± U.4o 


3 


no 1 n rr n 

9U ± U.5U 


2003 


Jun 22 




D.63 ± 0.57 


q trn _L not; 
o.oU ± U.zo 


o o o j_ n on 
o.oo ± U.zU 


4.57 ± U.7U 








2004 


Jan 31° 










/i 70 _l n on 


/I *7/1 _L n 0"7 

4. (A ± U.2/ 






2004 


oep 23 




A _1_ O Q Q 

4. DO ± 0.33 


O.U/ HZ u.zu 


0.4:0 Z±Z U.ZU 


AAA _L n OC 

4.44 ± U.zo 


/i 7n _l n /i n 

4. /y ± u.4y 





c:n _L n 40 


ZUuD 


Jan o 




4.4y ± U.oo 


2.91 ±0.16 


3.21 ±0.19 


/i m J_ n 07 


x ii in oq 
4. 11 ± U.zo 


3 


.18 ± 0.38 


900^ 


Apr 26 






2.58 ± 0.43 


2.92 ±0.17 


3 14 i 1 12 


Q 7Q _L n 00 






2005 


Jun 13° 






2.45 ±0.17 


3.11 ±0.19 


A ni _l n 04 


. < ^ u.oy 






2005 


Oct 10 












4.25 ± 0.55 






2006 


Mar 6 




3.99 ±0.26 


2.62 ±0.17 


3.20 ±0.16 


3.92 ±0.26 


3.93 ± 0.20 


2 


.91 ± 0.35 


2006 


Jun 6 






2.45 ±0.15 


3.20 ±0.17 


4.16 ±0.35 


3.61 ± 0.21 


2 


.32 ± 0.34 


2006 


Aug 30 e 






2.12 ±0.13 






3.24 ±0.43 






2006 


Dec 4 








3.13 ±0.31 










2006 


Dec 11 












4.54 ±0.91 






2007 


Apr 23 










3.69 ± 0.29 


3.57 ±0.20 


2 


.07 ±0.29 


2007 


Jul 31 






2.06 ±0.15 




3.77 ±0.32 








2007 


Aug 19 s 


8.9 ± 1.8 


3.07 ±0.13 




2.97 ±0.12 










2007 


Aug 23 h 












3.12 ± 0.19 


2 


.13 ± 0.23 


2008 


Oct 26 






2.27 ±0.12 













a Flux densities from fitting images and/or fitting u-v plane models. The tabulated uncertainties include a systematic 
uncertainty in the flux density calibration of 5% or more at frequencies of 22 GHz and below, and 10% or more at 43 GHz. 
They also include the estimated uncertainty of the galaxy subtraction for those epochs/frequencies at which the galaxy 
contribution was significant. We only tabulate measurements at resolutions < 7", and ones where the galaxy contribution 
was < 50% of the brightness of SN 1986J. 
b Observations between 4.86 and 4.99 GHz. 
c Observations between 8.41 and 8.46 GHz. 
d Observations between 22.21 and 22.66 GHz. 
c Data from the VLA archives. 

f Flux density at 1.7 GHz, rather than 1.43 GHz. 

s Values arc the average from two observing runs on 2007 Jul 31 and 2007 Sep 6. 
" Values are the average from two observing runs on 2007 Jul 31 and 2007 Sep 15. 



2.3. Radio Spectral Index 

As is clear from Figure [21 the radio spectral index, 
a is a strong function of both time and frequency. Up 
to 1995.5, however, the spectrum at least between 1.4 
and 15 GHz seems well described by a singl e power-law. 
Earlier on, up to 1988, iWeiler et al.l (|1990l ) reported a 
spectral index for the optically-thin part of — 0.671q'q4. 
IBall fc Kirkl ([l995t ). for the period up to 1991, find a 
somewhat flatter value of — 0.57 ± 0.01 from fitting a dif- 
fusive acceleration model. Fitting a power-law to the first 
set of flux densities included in Table [TJ namely those of 
1995.5, we find a notably flatter spectrum, with a spec- 
tral index of —0.44 ± 0.03 (we note that if we drop the 
22 GHz flux density from this fit on the suspicion that 
it might already be contaminated by the inversion, the 
fitted spectral index changes by less than the uncertain- 
ties). So the spectrum seems to have flattened signifi- 
cantly between 1988 and 1995. 

After 1995, the inversion appears in the spectrum, and 
its evolution becomes more complex. We will focus first 



on the low-frequency part o f the spectrum, which w e as- 
sociate with the shell. In [Bietcnho lz et al.l (|2002ft . we 
performed a two-part decomposition of the spectrum for 
the period 1998 to 2002, and obtained a spectral index 
of —0.55^_®'iQ. For even more recent times, we can cal- 
culate some representative spectral indices from the flux 
densities in Table Q] From our measurements on 2007 
Aug 19, we find = -0.71 ± 0.14. Although the 

uncertainties are rather large, this suggests also that the 
part of the spectrum at the lowest frequencies is again 
becoming steeper since 1995.5. 

In summary, the evidence suggests that the part of the 
spectrum below the inversion, which we associate with 
the shell, was fairly steep with a ~ —0.62 at times till 
1988, flattening to a ~ —0.44 in 1995, and once again 
steepening somewhat to a ~ —0.7 since then. By aver- 
aging the spectral indices derived by us and others for 
the low-frequency/optically-thin part of the spectrum, 
we find an average value over the whole timerange of 
—0.61 with an rms of 0.11. 

Another region of interest is the spectrum above the 
high-frequency turnover. This part of the spectrum is 
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Figure 1. Multi-frequency radio lightcurves for SN 1986J, as determined from VLA observations. We plot the lightcurves at 1.4, 5.0, 
8.4, 15, and 22 GHz. The flux-density scale on the left axis is relevant directly for the lightcurves at 1.5 and 5.0 GHz (marked with an 
asterisk), while those at 8.4, 15 and 22 GHz are shifted logarithmically downwards for better visibility, so that the shape of the lightcurve 
is preserved even though the flux density scale is not. The curves are shifted by the following factors: 8.4 GHz by 0.56, 15 GHz by 0.32, 
and 22 GHz by 0.14. The uncertainties are estimated standard errors, including statistical and systematic contributions. The data include 
both our own and re-reduced archival data (see Table [2}. 



dominated by emission from the central bright spot, and 
thus presumably represents its unabsorbed spectrum. A 
weighted average between 2006 Mar 6 and 2007 Aug 23 
gives aS = -0-61 ± 0.11. 

3. VLBI IMAGES 

In Figures [3] and U we show the VLBI images of 
SN 1986J at the different frequencies and dates in Ta- 
ble [U For comparison, we include in Figure [3] also some 
of our previously published images from 1999 and 2002. 
We discuss the images in turn, first those at 5 GHz, which 
together with the earlier images present the best picture 
of the evolution of the supernova, and then those at 8.4 
and 22 GHz. 

Since all the observations in this paper were phase- 



referenced to the brightness peak of 3C 66A, which is 
expected to be closely related to the core of the galaxy, 
we can accurately align the images of SN 1986J at differ- 
ent times and frequencies to facilitate inter-comparison. 
Unfortunately, since the earliest epochs of VLBI ob- 
servations of SN 1986 J (between 1987 and 1992, see 
iBietenholz et al.l 120021 ) were not phase- referenced, we 
cannot pinpoint the e xplosion location accur ately, as we 
did for SN 1993J (see IBietenholz et al.ll200l . However, 
we can visually estimate the position of the center of the 
shell sufficiently accurately for our purposes. We use this 
center position as the origin of our coordinate system in 
all the images, and also as our estimate of the explosion 
location. 

In particular, we estimated by eye the coordinates 
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Figure 2. The evolving radio spectrum of SN 1986 J, as determined from VLA observations. Each curve shows the radio spectrum at the 
epoch indicated at left, with the earliest spectrum at the top. The flux-density scale on the left axis is relevant directly for the first three 
curves from the top (1995.5, 1999.1 and 2002.4), whose epochs are marked with an asterisk. The other curves are shifted logarithmically 
progressively downwards for better visibility, so that the spectral shape is preserved even though the flux density scale is not. The curves 
are shifted by the following factors: 2004.7 by 0.80, 2005.0 by 0.60, 2005.3 by 0.50, 2006.2 by 0.25, 2006.4 by 0.15, and 2007.6 by 0.10. 
The uncertainties are estimated standard errors, including statistical and systematic contributions. The data include both our own and 
re-reduced archival data (see Table [2{. 



of the center of the shell in each of the four phase- 
referenced, 5-GHz images shown in Figure [3l We then 
averaged these four sets of coordinates and use the 
average to represent the coordinates of the shell cen- 
ter and the explosion center position, which is Ra = 
02 h 22 m 3K321457, decl. = 42° 19' 57/25951. We note 
that a minimum uncertainty in this center position can 
be derived from the scatter of the individual estimates, 
and is ~ 60 fias and ~ 180 /zas in RA and decl. re- 
spectively, and we estimate realistic uncertainties in the 
actual position of the center of the shell the order of 
200 /ias. We note that the largest component of this 
uncertainty is not astrometric, but comes from the diffi- 
culty of identifying the center position of the shell in the 



morphology, and our measurements of proper motion be- 
low do not depend on this estimate of the position of the 
shell's center. 

We discuss the images in turn: 

a) 5 GHz, 1999 Feb. 21: (Fig. Hi) The northeast bright 
spot dominates. There is a suggestion that it has an 
extension towards the cent er. A version of this image 
was previously published in iBietenholz et aT] ()2002[ ). 

b) 5 GHz, 2002 Nov. 10: (Fig. [3b) The northeast bright 
spot continues to dominate the image. Again, a slight 
extension towards the center is seen, but it cannot be 
clearly identified as the central bright s pot. A version of 
this image was previously published in IBietenholz et all 
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Figure 3. A sequence of VLBI images of SN 1986J at 5 GHz, showing its evolution over the last decade. In each panel, the 
coordinate system is centered on the estimated position of the center of the shell (02 h 22 m 31*321457, 42° 19' 57'.'25951, see 
text). The contours are in % of the peak brightness, with negative ones being dotted. Each panel has contours at 10, 20, 30, 
40, 50, 70, and 90%, with the 50% one being emphasized and the lowest white one. The lowest contours, given for each panel 
below, are three times the rms background brightness. The FWHM size of the convolving beam is indicated at lower left in 
each panel. North is up and east is to the left. This figure is also available as an mpeg animation in the electronic edition of the 
Astrophysical Journal, a) 1999 Feb 21, 5 GHz. The peak brightness was 1720 /j,3y beam -1 , the lowest contour 5.5%, the rms 
background brightness 31 /iJy beam - , and the convolving beam 2.09 x 1.02 mas at p. a. —7°. b) 2002 Nov 10, 5 GHz. The peak 
brightness was 543 /iJy beam -1 , the lowest contour 7.6%, the rms background brightness 14 /iJy beam -1 , and the convolving 
beam 1.57 x 1.12 mas at p. a. —10°. The visibility data were tapered slightly in u to increase signal-to-noise and produce a 
somewhat less elongated beam, c) 2005 Oct. 24, 5.0 GHz. The peak brightness was 414 /iJy beam -1 , the lowest contour 9.3%, 
the rms background brightness 13 /iJy beam -1 , and the convolving beam 1.65 x 0.83 mas at p. a. —20°. d) 2008 Oct. 25, 5.0 
GHz. The peak brightness was 388 /iJy beam -1 , the lowest contour 7%, the rms background brightness 9.0 /xJy beam -1 , and 
the convolving beam was 2.02 x 0.98 mas at p. a. —13°. 
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Figure 4. VLBI images of SN 1986J at 8.4 and 22 GHz. In each panel, the coordinate system is centered on the estimated position of 
the center of the shell (02 h 22 m 31?321457, 42° 19' 57'.'25951, sec text). The contours are in % of the peak brightness, with the lowest 
being at three times the rms background brightness, and the 50% being emphasized and being the first white contour. The FWHM of the 
convolving beam is shown at lower left. North is up and east to the left, a) Contours are drawn at —6, 6, 10, 20, 30, 50, 70, and 90% of the 
peak brightness, which was 717 ^Jy beam- 1 . The rms background brightness was 15 fijy beam 1 , and the convolving beam FWHM was 
1.24 X 0.68 mas at p.a. —12°. b) Contours are drawn at —17, 17, 30, 50, 70, and 90% of the peak brightness, which was 580 ^Jy beam -1 . 
The rms background brightness was 34 /zjy beam -1 , and the convolving beam FWHM was 0.85 X 0.35 mas at p.a. —20°. 



8 



Bietenholz, Bartel & Rupen 



(|2004aD . 

c) 5 GHz, 2005 Oct. 24: (Fig. [3fc) The central bright 
spot is now clearly seen at 5 GHz although it is not yet 
as strong as the northeast bright sp ot. A version of this 
image was previously published in iBietenholz fc Bartell 
(I2008an . 

d) 5 GHz, 2008 Oct. 25: (Fig. EH) The central bright 
spot now dominates at 5 GHz, and the northeast bright 
spot is becoming relatively fainter. The shell continues 
to expand (see § 13.11 below for a determination of the 
expansion velocity). 

e) 8.4 GHz, 2006 Dec. 3: (Fig. Ha) This run benefited 
from the newly available 512 Mbits s _1 recording band- 
width. The resulting image is dominated by the central 
bright spot, with the northeast bright spot and the re- 
mainder of the shell being only marginally discernible. 

f) 22 GHz, 2005 Apr. 25: This image was very simi- 
lar to the later 22 GHz image in 2006, discussed be- 
low. As it was of notably lower signal-to-noise ra- 
tio, we do not reproduce i t here (it was reproduced in 
IBietenholz fc Bartel|[2008al) . 

g) 22 GHz, 2006 Dec. 10: (Fig. Hb) At this frequency, 
the steep-spectrum shell is not visible, and the central 
bright spot dominates. The latter is now clearly resolved 
and appears somewhat asymmetrical, with a brightness 
peak to the east. The finger-like extension to the south 
is of marginal significance. This image was deconvolved 
with CLEAN, but we note that maximum entropy de- 
convolution (AIPS task VTESS) produces a very similar 
image. 

3.1. Expansion 

The average expansion velocity and deceleration can be 
determined from the VLBI data by estimating the size of 
the supernova at each epoch. As SN 1986J's morphology 
is becoming increasingly irregular, estimating the size by 
fitting a simple geometrical mod el to the visibility data 
(as we did fo r SN 1993J, see e.g.. IBietenholz et al.l[200ll : 
Bartel_elaL|[2002) is no longer practical. We take as a 
representative (angular) radius the average radius of the 
contour which contains 90% of the total flux density in 
the image. Measures of the radius derived from images 
are complicated by the convolution of the images by a 
clean beam, but a consistent estimate of the expansion 
rate can be obtained by comparing radii derived from im- 
ages convolved with an app roximately co-evolvin g beam. 
For a fuller discussion, see IBietenholz et al.l ()2002f ) . where 
we showed that the expansion up to 1999 had a power- 
law form with the radius, r oc i°- 71 ± 011 . We accordingly 
use images convolved with a co-moving beam0 whose di- 
mensions scale with t - 70 . We call this radius # 90 % nux - 

6 This is the same procedure we used inlBTetcnholz et al. (2002), 
where we also found that the expansion curve is not sensitive to 
the exact time-dependence of the restoring beam, as long as it is 
similar to the time-dependence of the radius. The use of an approx- 
imate value for the coefficient of the expansion power-law should 
therefore not significantly bias our result. Note that the value of 
0.70 we use for this co efficient is slightly diffe rent from the value 
of 0.75 that we used in Bietenholz et al. (2002). The present value 
is a somewhat better match to the fit power-law expansion of the 
supernova. Note also that our present value of the expansion coef- 
ficient is in agr eement within the uncer tainties with the value we 
had obtained in Bietenholz et al. (2002] ), but does not support th e 
higher value of 0.90 ± 0.06 reported by Perez- Torres et al. (2002). 



Table [3] gives the derived radii for our various observing 
epochs. For observations after 2000, we do not determine 
radii from the observations at 8 GHz or above because 
at these frequencies, the shell was not sufficiently dis- 
tinct for a useful measurement. We further note that the 
increasingly non-self-similar evolution of the morphology 
makes such average measures of the radius and expansion 
less reliable, however, our fit to the expansion should give 
a reasonable estimate of the properties of the expansion 
averaged over the entire shell. 

We fix the explosi on date, to, at 1983.2, the value 
that we determined in IBietenholz et all ([20021 ) from both 
the expansion measurements and the radio light-curve. 
A least-squares fit of a power-law then gives an av- 
erage expansion 9O % a ux oc £0- 69 ±°o 3 ; consistent with, 
but more accurate th an the value we had derived in 
IBietenholz et al.l (|2002t ) for the measurements up to 1999. 

We plot the resulting expansion curve in Figure [5j The 
two radii after epoch 1999 fall slightly below the earlier 
power-law fit, suggesting a possible slight increase in the 
deceleration. However, we caution that the increase in 
brightness of the central bright spot relative to that of 
the shell might lead to a #90% flux being a slight underes- 
timate of the true radius of the shell. 

The average expansion velocity measured between the 
1999 and 2008 epochs was 5700±1000km s" 1 . Although 
this velocity, based on 9qq% a ux , is likely not the exact 
velocity of the outer radius of the shell of radio emission, 
it should be a reasonable approximation thereof. 

3.2. Positions and Proper Motion of the Two Bright 

Spots 

How does the expansion of the supernova compare with 
the proper motions of the northeast and central bright 
spots? Since our images are phase-referenced to the 
brightness peak of the un-related nearby compact radio 
source, 3C 66A, we can accurately determine the proper 
motion of the two bright spots relative to 3C 66A. 

Before we proceed to the proper motions, we briefly 
discuss the uncertainties which will apply to such mea- 
surements (for other discussions on the astrometric ac- 
curacy obtainable with VLBI, see IBietenholz eTaIll200l 
l2004bl : IPradel et al.l 120061 ). Our estimated astrometric 
uncertainties below consist of the following three compo- 
nents: 1) An uncertainty due to the noise in the image. 
2) An uncertainty in identifying a particular point in the 
morphology of the reference source, 3C 66A, which is 
slightly time and frequency dependent^. We estimate the 
latter uncertainty component as 20 /zas in RA and 70 /ias 
in deck, with the deck component being larger because 
the source is elongated approximately in the north-south 
direction. 3) An astrometric uncertainty due to errors 
in the station coordinates, earth orientation, and tropo- 
sphcric correcti on, which we estim ate at 30 //as in each 
coordinate (see IPradel et al.l[2006l ). 

We first turn to the northeast bright spot. It was 
clearly detected already in 1999 (Bietenholz et al., 2002; 
for a similar image from an independent reduction of 

7 IKovalev et al.l 1 120081) show that frequency-dependent shifts in 
core positi on are common among co mpac t radio sources. For ex - 
ample, see IBietenholz et al.l l|2004bl 1 and [Bietenholz et al. ( 2000) 
for measurements of the variability in time and frequency of the 
compact core of the nearby galaxy M81. 
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Table 3 

Angular Sizes of SN 1986J 



Date 


Frequency 


Age a 


Angular Size b 




(GHZ) 


(yrs) 


(mas) 


1988 Sep 29 


8.4 


5.55 


1.37±0.10 c 


1990 Jul 21 


8.4 


7.35 


1.83±0.10 c 


1999 Feb 22 


5.0 


15.94 


3.07±0.10 c 


2002 Nov 10 


5.0 


19.66 


3.60 ±0.14 


2005 Oct 24 


5.0 


22.62 


3.71 ±0.15 


2008 Oct 26 


5.0 


25.62 


4.23 ±0.17 



a The explosio n date is not pre c isely known. We use a best-fit explosion epoch of 1983.2 
determined in Bictcnhofz ct al. ( 20021) . 

b We take as a representative angular radius 9 S0 % Hux , which is the average radius of the 
area containing 90% of the clean flux density in the image. The images are convolved 
with an approximately co-moving restorin g beam. See text f or de tails. 
c These values are from data discussed in iBietenholz et"aTI (|2002l ). although the present 
values of 9 aQ % a ux differ slightly from the ones reported there because we are using a slightly 
different co-moving convolving beam. In that earlier reference, we used a convolving beam 
whose size increased oc t 0,75 , whereas the present values are based on a beam whose size 
increases oc t 0,70 , which better matches the actual expansion of the supernovae. The 
difference in derived radii are small, with the values of g0 % flux changing by less than 
their stated uncertainties. 
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Figure 5. The expansion curve of SN 1 986J . We plot the angular outer radius, t9go% fl ux (the angular radius which contains 90% of the 
total flux density in the image, see text, § I3.1l for details), against the age, both on logarithmic scales. The line indicates expansion of the 
form 8 oc (i — io) m , with m = 0.69, where t = corresponds to the fixed epoch 1983.2, which was the best weighted least-squares fit to the 
points plotted. 



our observations, see IPerez-Torres et al.l l2002). We find 
that, based on a linear fit to the position in the 1999 
Feb., 2002 Nov., 2005 Oct., and 2008 Oct. images at 
5 GHz, the proper motion of the northeast bright spot 
was 60 ± 20 /ias yr^ 1 , corresponding to a speed of 
3000 ±1000 km s"\ 2400 ±250 and 1500 ± 700 km s" 1 , 
in RA and decl. respectively. 
Does the current proper motion of the northeast bright 



spot place it at the explosion center at the time of ex- 
plosion? To answer this question requires knowing the 
position of the explosion enter. As mentioned earlier, 
this center is not accurately known for SN 1986J, but 
as an estimate we use the average position of the cen- 
ter of the shell, as described in §[3] above. Extrapolating 
the northeast bright spot's present proper motion the ex- 
plosion epoch of 1983.2 yields a position consistent with 
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the estimated center position within 1.9 and 0.8 sigma in 
RA and decl., respectively. If allowance for deceleration 
is made, then the spot's 1983.2 extrapolated position is 
even closer to the center. In fact, a radial motion with 
the same power-law dependence on time as the overall 
expansion, which has r oc i 069 , is compatible with our 
measurements, as can be seen in Figure[5J where we show 
the radial motion of the northeast bright spot. 

A homologous, power-law expansion which includes 
the northeast bright spot is therefore consistent with 
our measurements, with the spot's present position and 
proper motion being compatible with it being at the 
center of the shell at the explosion epoch. We note, 
however, that the uncertainties are not small, so consid- 
erable deviation from the power-law expansion, for ex- 
ample stronger deceleration, or motion with a constant 
speed, are also consistent with our measurements. 

We turn now to the central bright spot. It is best 
defined in the images from 2003 June 22 (15 GHz, 
iBietenholz etafl2004a| ). 2006 Dec. 10 (22 GHz) and 2008 
Oct. 26 (5 GHz). If we estimate the center bright spot's 
position at each epoch by the position of the brightness 
peak, we can determine the spot's proper motion. A lin- 
ear fit to the three positions gives a proper motion of 
+4 and +31 /j,as yr -1 in RA and decl., respectively. The 
morphology of the bright spot makes the peak brightness 
position somewhat resolution-dependent, (and possibly 
frequency-dependent), and we estimate that the accu- 
racy with which any particular spot in the morphology 
can reliably be identified to be on the order of 100 /ias. 
In addition to the other uncertainties estimated above, 
we use a rounded value of 120 ^as for the uncertainty in 
the spot position in each coordinate, which leads to an 
uncertainty in the proper motion of 31 /ias yr -1 in each 
coordinate. This implies a speed of 1500±1500 km s _1 at 
10 Mpc. The estimate of the proper motion is consistent 
with the central bright spot being stationary. The aver- 
age position of the central bright spot is 300 ±150 /xas 
from the geometric center of the shell (see § [3]). How- 
ever, given the uncertainty in associating the explosion 
center with the geometric center of the shell, we only say 
that the position of the central bright spot is marginally 
consistent with that of the explosion center. 

3.3. Brightness Evolution of the Two Bright Spots 

In order to quantify the evolution of the brightness of 
the northeast and central bright spots, we fitted simple 
models to each bright spot in the image plane. In par- 
ticular, in the region near each spot we fit an elliptical 
Gaussian of the same dimensions as the restoring beam, 
i.e., the equivalent of a point source, and a baseline level. 

Although the assumptions that the bright spots are in 
fact unresolved and that they can be separated in this 
way are not unquestionable, we feel that this procedure 
gives a reasonable approximation of the relative flux den- 
sity of the two spots at each epoch, and is more accurate 
than, for example, relying only the spot's peak bright- 
ness. 

In Figure [7] we show the 5-GHz evolution of the flux 
density of each bright spot relative to the total flux den- 
sity in the image as a function of time. The north- 
east bright spot's relative brightness increases, reach- 
ing a maximum of >20% of the total flux density at 
an age of ~16 yr, and fades thereafter. Recall, how- 



ever, that the total flux density is decreasing relatively 
rapidly as a function of time (§ 12.21 Figure [lj, so in 
absolute terms, the brightness of the northeast bright 
spot decreases monotonically. At its relative peak, its 
flux density was ~ 1.3 mJy, and its radio luminosity, 
assuming an a = —0.61 spectrum up to 100 GHz, was 
- 6 x 10 36 erg s _1 . 

The central bright spot was not clearly visible (at 
5 GHz) till age 22.6 yr, but by age 25.6 yr, it had in- 
creased to ~15% of the total flux density. The central 
bright spot is still increasing in brightness in absolute 
terms, but the absolute increase is not as rapid as the 
relative one plotted in Figure [7] In 2008, the central 
bright spot had a flux density of ~ 390 fxJy, and a ra- 
dio luminosity of ~ 2 x 10 36 erg s _1 (again assuming a 
spectrum with a = —0.61 up to 100 GHz). 

4. DISCUSSION 

4.1. Overall Evolution: Expansion Curve, Deceleration, 
and Spectral Index 

The expansion curve in Figure [5] is compatible with 
a power-law expansion, and earlier we found the best 
fit to correspond to r oc t o eg . The flux density decay, 
however, shows a distinct change in the power-law slope 
at some time between 1989 and 1995 , i.e., between ages 
of 6 to 12 yr (Bietenho lz et alj |2002). Unfortunately no 
flux density measurements are available in this period 
to determine the time of change more accurately. Simi- 
lar behavior was found for SN 1993 J with a steepening 
of the flux d ensity decay at i ^7 yr (|Bartel et alj|2002t 
iWeiler et al1l2007l) . In the case of SN 1993 J, this steep- 
ening was found to corre spond to a decrease in deceler- 
ation (Bartel et al. 2002). This is the general behavior 
expected when there is a relative decrease in the CSM 
density. Is a similar behavior seen in SN 1986J? Exam- 
ining the expansion curve in Figure [5] (see also expansion 
curve s including earlier points in IBietenholz et al]|2002L 
2005) suggests that a significant decrease in deceleration 
in the age range of 6 to 12 yr is not compatible with the 
data, although a small increase in deceleration might 
be. This would suggest that the increase in the rate of 
flux density decay is due not to a relative decrease in 
the CSM density but rather to one in the ejecta density. 
This suggests a steep density profile in the ejecta in the 
exterior regions, with a somewhat flatter profile in the 
interior. Such a flattening is indeed expected f rom mod- 
els of supergiant stars (see, e.g., discussion in iChevalierl 
l2005f h Although suggestive, we do not regard this result 
as definitive because the measurements do not demand 
any break in the power-law expansion curve. 

The evolution of the integrated spectrum is complex, 
but at low frequencies the spectrum seems to have a 
power-law shape. One can with reasonable confidence 
associate the part of the spectrum at the lowest frequen- 
cies with the shell, with the complex evolution due to 
the emergence of the central bright spot being restricted 
to higher frequencies (> 5 GHz). As we showed in § 12.31 
above, the lowest part of the spectrum seems be initially 
rather steep, with a = —0.67, then by 1995 it flattens 
somewhat to a = —0.44, and subsequently it steepens 
again. 

What could cause such a variation of the radio emis- 
sion spectrum with time? Since the emission spectrum of 
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Figure 6. The angular distance of the northeast bright spot from the center as a function of time. The additional axis labels along the top 
give the cal enda r year. The spot's angular distance is taken to be that between the spot local maximum and the estimated explosion c ente r 
(see text, § 13.21 for details on the estimation of the explosion center). The shell expands with outer angular radius oc t - 69 ( see text § 13.111 . 
Of the possible expansion curves which have 9 oc t - 69 , we show as the solid line the one that best matches the spot's motion. Note that the 
spot's displacement as a function of time is somewhat dependent upon the choice of the explosion center, although for reasonable choices, 
the nature of the expansion curve remains unchanged. Our plotted error bars are the uncertainty in the spot position and an estimated 
uncertainty of 0.1 mas in the explosion location added in quadrature. The spot's motion is consistent with homologous expansion. 
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Figure 7. The evolution at 5 GHz of the fractional flux density of the bright spots in the northeast and in the center of the projected 
shell as a function of time. The additional axis labels along the top give the calendar year. The bright spot flux densities are estimated by 
fitting, by least squares, a point source and a baseline level in the region near the spot. The fraction of the flux density is taken to be the 
fitted flux density for the point source divided by the total CLEAN flux density in the image. For the first three epochs, the bright spot 
in the center is not discernible. The chevron shows an upper limit on the central bright spot flux density at age 19.7 yr, which we derived 
by taking half the surface brightness at its location. 
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synchrotron radiation reflects the energy spectrum of the 
relativistic particles, one obvious interpretation is that 
the latter energy spectrum is also changing with time, 
perhaps in response to varying conditions at the shock 
front where the particle acceleration occurs. However, 
optical depth effects also offer a possible explanation. Al- 
though most of the radio-emitting shell is optically thin 
at cm wavelengths, the interior of the shell is not ex- 
pected to be so because of the high density of thermal 
material. If mixing is occurring, a substantial transition 
zone between the optically thin shocked material and the 
optically thick ejecta may exist, which could partly ab- 
sorb the radio emission from a substantial fraction of the 
shell, and thus modify the integrated spectrum. In fact, 
the slow turn-on of the lightcurve suggests a significant 
amount of absorption mixed in with the radio emission 
already at early times (|Weiler et aLlll990h . 

We note that changes with time in the radio spec- 
tral index are seen in other supernovae. In particular, 
a flattenin g with time of the radio spectrum as seen for 
SN 1993 J ([Bartel et al.l^OOafl and at least at the lat - 
cst times also for SN 1979C (jBartel fc BietenhollllM h 
may be characteristic of the evolution of radio super- 
novae, and could be coupled with a steepening of the 
lightcurve decay. 

4.2. Comparison of Expansion Curve to Optical Radial 

Velocities 

From our VLBI measurements, we determined the ex- 
pansion velocity of the edge of the radio emission region 
in the plane of the sky. Since the forward shock forms 
the outer boundary of the region where radio emission 
is expected to occur, our VLBI expansion velocity can 
be reasonably associated with the velocity of the for- 
ward shock velocity (see discu ssion on this subject for 
SN 1993 J in lBartel et~al~ll2007t ). Note that since the ge- 
ometry of SN 1986J is somewhat irregular, the assump- 
tion of sphericity which allows a comparison of the ve- 
locities in the sky-plane measured with VLBI with those 
along the line-of-sight measured from optical spectra is 
likely true only to first order. 

In the optical, there are a number of spectral lines 
present which allow radia l veloc ity measurements. Re- 
cently, iMilisavljevic et all ((2008), obtained new optical 
spectra at epoch 2007.7, and compared them to earlier 
ones from 1989.7 and 1991.8. The optical emission lines 
are thought to mainly come from two distinct locations. 
Narrow optical lines, chiefly Ha, He I and [N II] A5755, 
are thought to originate in dense, shock-heated circum- 
stellar medium, and therefore have relatively narrow line- 
widths, typically < 2000 km s" 1 . These lines were more 
prominent early on and have faded since 1989. 

Currently the most prominent emission lines are the 
[O I] AA6300,6364 and [O II] AA 7319, 7330 forbid- 
den lines of oxygen, with the [O III] AA4959, 5007 line 
also being discernible. In 1989.7, these lines showed a 
double-peaked structure, with one peak at a velocity of 
— 1000 km s _1 (relative to the systemic one), and the 

8 ICha ndra ct al. ( 2004]) also find a flattening with time of the 
optically-thin spectrum of SN 1993J. A slight but systematic 
flattening is also v isible in the plotted data after t = 500 d of 
IWeiler ct al. (2007), although they fitted a model with a constant 
spectral index for the optically thin part of the spectrum. 



other near —3500 km s , and the maximum velocity, on 
the blue sides of these lines, was ~ —6000 km s _1 . Simi- 
lar velocities were seen in 1991.8. By 2007.7, the double 
structure was less clearly discernible, the maximum blue 
velocity h ad decreased and the [O III] line had almost 
vanished. Milisavlievic et al . (2008) attribute these for- 
bidden lines to ejecta which has been heated (but not yet 
shocked) by the reverse shock. 

This interpretation of the optical lines is consistent 
with the expansion velocities measured in the radio. In 
1989.7, the forward shock velocity was ^8300 km s _1 , 
and the reverse shock velocity would be expected to be 
about 80% of this, or ~6600 km s _1 . The velocities im- 
plied by both the line-center and the total widths seen 
in the forbidden oxygen lines are consistent with being 
inside the reverse shock in 1989.7. By 2007.7, the for- 
ward shock velocity is ^5500 km s _1 , and the expected 
reverse shock velocity ~4400 km s _1 . The evolution of 
the optical line profiles is consistent with this reduction 
in the velocity of the reverse shock in that, particularly 
for the [O II] line, with the highest (blue) velocities now 
being ~ —4300 km s _1 . This suggests that, of the ejecta 
responsible for the oxygen forbidden line emission, the 
highest velocity fraction passed through the reverse shock 
between 1989.7 and 2007.7, and the remainder is now sit- 
uated relatively close to the reverse shock. 

4.3. Magnetic Field 

The magnetic field can be estimated from the radio 
brightnes s and the source si ze using minimum energy ar- 
guments (jPacholczvkl Il970f1. Given that the spectrum 
shows deviations from a power-law, we will perform the 
calculation only for the shell component, which domi- 
nates below v = 5 GHz and is observed to have a power- 
law spectrum with an average a = —0.5 (§ I2.3[) . The 
derived field strength depends on the ratio, fc, of ener- 
gies of relativistic protons and electrons. We use a value 
of fc = 100, as was found to b e reasonable for stron g 
shocks in supernova remnants bv lBeck fc KrausH (|2005|) . 
We estimate first the field strength in the shell, for which 
we use the angular radii derived in § 13.11 above, and we 
take the flux density at 1.4 GHz to be representative of 
the shell (as opposed to that of the central bright spot). 

Using the values from 1988 to 2008, we find that B 
has an approximate power-law dependence on time and 
radius, as might be expected from the approximately 
power-law evolution of both the flux density and radius. 
In particular, we find that for minimum energy, B ~ 
60 (t/lOyr)" 13 mG, and that B cx r" 18 for t > 1500 d. 
Although such estimates of the field are uncertain by fac- 
tors of several, the estimates of its dependence on t and 
r are likely to be somewhat better constrained with the 
exponents being uncertain by perhaps ^70%. 

4.4. The Mass Loss Rate of the Progenitor 

9 We note that Beck & Krausc ( 2005) point out a number of 
shortcomings in the standard formulae of Pacholczyk and give re- 
vised formulae. However, for relatively flat spectra, the differences 
between the standard and revised formulae are not large and only 
weakly dependent on the field strength, so a field calculated using 
the revised formula would have almost the same dependence on 
time and radius. In the interests of easier comparison with earlier 
resul ts we estimate our magnetic field using the standard formula 
from lPacholczvkl 119713) . 
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Can we estimate the mass-loss rate of the progeni- 
tor from our observations? This mass-loss rate is often 
estimated from studies of the radio lightcurves. It is 
determined mostly by the rising part of the lightcurve, 
while the emission is optically thick, and is usually ex- 
pressed as M/v w , where v w is the wind velocity. How- 
ever, such calculations can significantly overestimate the 
value of M/v w (see, e.g., iBartel fc Bietenhold [2003h . 
and in general have been shown to be somewhat prob- 
lematic if not based on com plete physical models (e.g., 
iFransson fc Bjornssonl 12005). In particular, in the case 

of SN 1986 J. iWeiler et all (|1990T ) cited a value of M/v w 
= 2.4xl0" 4 M Q yr _1 (for u w = 10 km s" 1 ), based not 
on the rising part of the lightcurve, but rather on a short 
term variation in the op tical depth seen in late 1988. 
Later. IWeiler et all (|2002f) reported a notably lower value 
of 4.3xl0 -5 Mq yr" 1 , which is base d on a model i nvolv - 
ing clumpy absorption in the CSM. iHouck et al.l (|1998f ) 
suggest a value of a few times 10~ 4 Mq yr -1 based on X- 
ray observations. In any case, given the departures from 
self-similarity and the considerable evidence of dumpi- 
ness in the CSM, which will have a strong effect on both 
radio brightness and absorption, any estimated mass-loss 
rate is likely to be quite uncertain. 

A mass-loss rate between 4xl0 -5 and, say, 3xl0" 4 
Mq yr -1 with v w = 10 km s _1 would imply that by age 
25.6 yr (2008.8), when the radius of the shock- front is 
6.3xl0 17 cm, the swept-up mass would be 0.8 ~ 6 Mq. 
For this range in mass-loss rate, the kinetic energy of 
the swept-up mass would be (0.2 ~ 1.7) x 10 51 erg 
(v ~ 5400 km s" 1 ). Since the kinetic energy of the swept- 
up mass can be at most a fraction of the original kinetic 
energy of the ejecta, which was likely not much larger 
than 10 51 erg (e.g., Uones et"al"1 I1998D . we regard the 
higher end of this mass-loss range as unlikely on energetic 
grounds. We therefore estimate a present swept-up mass 
of 0.8 ~ 2 Mq, and an average mass- loss rate (M/v w ) in 
the range of (4 ~ 10) x 10" 5 Mq yr" 1 . Such a range of 
mass-loss rate implies a CSM density at the 2008.8 outer 
shock radius of 6.3xl0 17 cm of (5 ~ 13)xl0~ 22 g cm" 3 . 
If we assume that the fraction of hydrogen atoms is 75% 
of the total, then the hydrogen number density is be- 
tween 230 ~ 570 cm -3 . 

4.5. The Northeast Bright Spot 

The images are characterized by two bright spots, 
whose flux densities vary. The first one, in the north- 
east, appears to be a component of the shell. It was first 
clearly discernible at age ~16 yr, although it may have 
been present already earlier. Its proper motion is con- 
sistent with homologous expansion with the bulk of the 
shell (Fig. [6|). As can be seen in Figure [3 at its peak 
relative brightness (age ~16 yr) approximately ~20% of 
the total flux density originates from this spot. Unfor- 
tunately, our observations do not constrain the relative 
lightcurve around the time of the peak precisely, the peak 
relative brightness may have been as high as ~30% of 
the total, although it must have occurred between our 
1990 and 1999 epochs. The northeast bright spot ap- 
pears marginally resolved in our images. On 2005 Oct. 24 
our FWHM resolution was 1.65 x 0.83 mas, so the spot's 
angular diameter must have been ~1 mas, corresponding 



to 14% of the shell's outer diameter, or 1.5 xlO 17 cm. 

Since the northeast bright spot appears to be part of 
the projected shell, the most plausible interpretation is 
that it is due to the ejecta interacting with a dense clump 
in the CSM. This would account for the relative bright- 
ness of the spot first increasing and then falling off. In 
fact, the relative lightcurve of the northeast bright spot 
in Figure[7]suggests that the shock has traversed the bulk 
of the clump between ages 10 to 20 yr. Our expansion 
curve suggests that, on average around the circumfer- 
ence, the shock travelled a distance of 2xl0 17 cm during 
this time, in reasonable agreement with the clump diam- 
eter estimated above. 

If we assume both the dense clump and the outer shock 
to be spherical, with the clump having the dimensions es- 
timated above, then, when the shock is midway through 
the clump, ~1% of the outer shock front area is inter- 
acting with the dense clump. Since the radio emission 
is driven by the shock, we will consider the radio lu- 
minosity per unit area of the shock, and we will further 
assume that the radio luminosity scales as the 5- GHz flux 
density. The northeast bright spot at its peak emitted 
about one quarter of the 5- GHz flux density of the re- 
mainder of the shell, we can therefore calculate that the 
region in the dense clump must be ~25 times brighter 
than the remainder. Our proper motion measurements 
show that the velocity of the shock in the clump is 
approximately the same as it is for the remainder of 
the shell, so the radio luminosity per unit area of the 
shock should scale with the CSM densitjB We found 
above f§ I4.4[) that the present average CSM density was 
(5 ~ 13)xl0 -22 g cm -3 , suggesting that the densities in 
the clump are (1.3 ~ 3.2) xlO -20 g cm" 3 , corresponding 
to hydrogen number densities of 5800 ~ 14,000 cm" 3 , 
giving the clump a total mass of 0.01 ~ 0.03 Mq. We 
note, however, that these densities represent an average 
over the clump, and much higher densities in parts of the 
clump are not excluded by the radio data. 

There are a number of other indicators pointing to 
the existence of such dense clumps in the CSM of 
SN 1986 J. First , the very high Balmer decrement 
(iRupen et al.l I1987T) . and other fe atures in the late-time 
spectra (jMilisavljevic et all 120081 ). suggest the presence 



of very high number densities (> 10 6 cm" 3 ). Sec- 
ond, the slow rise of the radio lightcu rve is inter- 
prete d as also implying a very clumpy CSM (jWeiler et al.l 
1990). Third, because the Ha line, which has a width 
of only ~10QQ km s" 1 (e.g., iMilisavlievic et al.l 120081: 
iLeibundgut et al.l fl991). is far too narrow to arise from 
the ejecta near the forward shock (v ~ 5700 km s ; 
§ 13. ip . it has been suggested that most of the Ha emission 
is coming not from the ejecta, but rather from shocked, 
dense clouds in the CSM resu lting from a very clumpy 
wind from the progenitor star fChugai fc Danzigerlll994l : 

10 The radio lumi nosity is though t to sc ale with the post-shock 
energy density (e.g., Chevalier 1998, 1982). In the absence of ab- 
sorption, the radio luminosity will depend on the fractions of the 
post-shock energy which go into the magnetic field and into accel- 
erated particles. There is no solid ground for thinking that these 
fractions are constant, indeed they are likely to vary somewhat 
both as a function of shock velocity and the density and chemical 
composition of the CSM. It is likely, however, that these variations 
are relatively small, whereas the variation in the CSM density can 
be several orders of magnitude or more. 
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IChugailH993h . 

Red supergiant stars of mass 20-40 Mq are known to 
have episodes of strong mass loss, resulting in large den- 
sity contras ts in the CSM. For example, iSmith et al.1 
([20011120091) show that the red supergiant VY Canis Ma- 
joris (mass 20-40 Mq), which is likely similar to the pro- 
genitor of SN 1986J, is surrounded by a nebula produced 
by rapid mass loss (2 ~ 4 x 10 -4 Mq yr _1 ). This neb- 
ula, which extends out to a radius of > 10 17 cm, similar 
in scale to the CSM required for SN 1986J, is highly 
structured with knots and filament ary arcs, wh i ch are 
distributed in an asymmetric way. ISmith et al.l (2001) 
argue that the structures in the wind were produced by 
ejection events which were localized on the stellar sur- 
face. For VY Canis Majoris then, and by implication 
also for other red supergiants, the evidence seems to sug- 
gest that the winds can be highly structured, with large 
density contrasts on small spatial scales. Exactly such 
structure seems to be required to reproduce the radio 
emission from SN 1986J. 

4.6. The Central Bright Spot: Is It Emission from the 
Compact Remnant of the Supernova Explosion? 

What is the nature of the central bright spot — is it 
the emission associated with a neutron-star or black-hole 
compact remnant of the supernova explosion? The first 
question would be whether the central bright spot is po- 
sitionally coincident with the explosion center. As we 
showed in § 13.21 above, the position of the central bright 
spot is in fact compatible with being in the center of 
the shell, and its proper motion is consistent with be- 
ing stationary with respect to the shell. Note that the 
expansion speed of the shell is much larger than the pe- 
culiar velocities of stars and even of pulsars, which latter 
have a veloci ty dispersion of 200 ~ 300 km s -1 in the 
Gala xy (e.g., iHobbs et al.l 120051 iLvne fc Graham- Smith! 
1990j), so the effect of any peculiar motion of the pro- 
genitor or any "kick" received in the explosion would be 
small and the compact remnant would be expected to re- 
main essentially in the center of the shell over the ~25 yrs 
since the supernova explosion. In this respect, our posi- 
tion and proper motion measurements for the central hot 
spot are consistent with the expectations for a black-hole 
or neutron star compact remnant. 

We first noted the inversion in t he integrated radio 
spectrum in IBietenholz et al.l ((2002), and we showed in 
iBietenholz et al.l (|2004ah that the inversion in the inte- 
grated spectrum was in fact associated with the central 
bright spot, with the shell and the northeast bright spot 
discussed above having a simple power-law spectrum as 
is usually seen in supernovae. In § 12.21 and Figure [2] 
above we show that both the inflection and the high- 
frequency turnover points in the radio spectrum are pro- 
gressing downwards in frequency as SN 1986J ages. The 
inverted spectrum of the central bright spot suggests 
partly absorbed radio emission. An opacity decreasing 
due to the expansion would cause the observed progres- 
sion to lower frequencies of both the inversion and high- 
frequency turn over points. 

We argued in IBietenholz et afl (|2005| ) that the absorp- 
tion is probably due to free-free absorption rather than 
to synchrotron self-absorption (SSA). In fact, our latest 
22 GHz images, where the central bright spot is some- 
what resolved, largely rules out SSA on the following 



grounds. Examination of the image (Fig. 0t>) suggests a 
size of > 1 x 10 17 cm for the central bright spot. The 
radius below which SSA would be im portant, R p , can be 
calculated following iChevalierl (|1998l ): 



R p = 



8.8 x 10 15 k 



-1/10 f / 
0.5 



-1/19 



9/19 



VM PC y 



18/19 



(sGHz) 



where F p is the flux density at the spectral peak, k is the 
ratio of relativistic electron energy density to magnetic 
energy density, D is the distance, and / is the filling 
factor. Assuming equipartition (k — 1) and substituting 
our values of F p = 3.8 mJy at u p = 14 GHz (2007.6), 
we obtain Rp ~ 1 x 10 16 cm, which is far lower than the 
size determined from the image of > 1 x 10 17 cm, thus 
effectively ruling out significant SSA at our observing 
frequencies. 

The appearance of the central bright spot at a time 
long after the supernova's integrated radio spectrum was 
optically thin at cm wavelengths suggested that any re- 
maining absorption due to the circumstellar material was 
minimal. It was natural, therefore, to assume that the 
absorption seen for the central bright spot was due to 
un-shocked material interior to the shocked shell, rather 
than to the CSM. With this assumption, the central 
bright spot's likeliest physical location was in the cen- 
ter of the shell, and its most likely interpretation was 
in terms of radio emission associated with the neu- 
tron star or black hole e xpected to have been left be- 
hind after the explosion ( [ Bietenholz et ai]l2004al 120051 : 
IBietenholz fc Bartelll2008a|) . 

The turnover frequency is approximately equal to the 
frequency at which the optical depth to free-free absorp- 
tion, r, is unity, is given by 

i/ T =i = 0.3 {T-^Nl dl) 1 ' 2 GHz, 

where T e is the electron temperature in K, N e is the num- 
ber density of electrons in cm" 3 , assumed constant along 
the path-length in pc, dl, with N 2 dl being the emission 
measure. Our integrated spectra (Fig. [5]) show that the 
spectral peak due to the central bright spot occurs at 
~20 GHz in 2002.4. If we assumeEJ a T e of 10 4 K then 
Nidi is 1.1 xlO 9 cm" 6 pc. By 2007.6, the turnover fre- 
quency is ~14 GHz, so Ag dl = 0.6 xlO 9 cm -6 pc. 

The temporal dependence of the absorption is con- 
sistent with that expected from the gas within the ex- 
panding shell. We found that the emission measure de- 
creased by a factor of ~2 between 2002.4 and 2007.6. If 
we assume the system is expanding homologously, then 
the density within the shell is proportional to r~ 3 and 
Ag dl oc r~ 5 . If we further take r oc t a e9 as we found for 
the shell, we would expect dl to be oc i -3 - 45 , which 
would lead to a decrease by close to the observed factor 
of 2 between 2002.4 and 2007.6. 

We note, however, that the inverted part of the spec- 
trum is much flatter than is expected from free-free ab- 
sorption, which produces spectra with an exponential 

Although T e is not well known, and can vary strongly with 
radius, th is is probably a reasonable v alue for radii > 1 X 10 17 cm 
(see, e.g., Lundqvist & Fransson 1988). 
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cutoff at low frequencies. Such flatter spectra below the 
turnover are produced when a range of different optical 
depths is present. This suggests that the absorbing ma- 
terial is somewhat fragmented, with some lines of sight 
having much lower optical depths than others. Indeed, 
fragmentation of the ejecta is not surprising, since there 
are a number of insta bilities ope r ating in the expand- 
ing shell of eiecta fe g.. lGulll ll973H Jun fc Norman! 1 19961: 
IChevalier fc Blondinl ll995 N ) , which are expected to lead 
to fragmentation. If the central bright spot is due to a 
pulsar, then further instabilities occur when the young 
pulsa r's wind nebula expands into the supernova eject a 
(e.g. JBandiera et alJ[l983HChevalier fc Franssonlll992ft . 

As we have shown, the absorption of the radio emis- 
sion from the central bright spot is consistent with what 
might be expected from the intervening material in the 
shell and its observed expansion rate. This suggests that 
the integrated spectrum above the turnover point shows 
the intrinsic, un-absorbed spectrum of the central bright 
spot. In fact, the integrated spectra in Figure [2] show 
that the spectral index above the turnover point is very 
similar to that below the inversion point. So, the intrin- 
sic spectral index of the central bright spot is also similar 
to that of the shell. 

If the central bright spot were in fact in the physical 
center of the shell and represented radio emission asso- 
ciated with the compact remnant of the supernova - 
either a black hole or a neutron star — then it would be 
somewhat of a coincidence for its radio emission to have 
a similar spectral index as that of the shell emission. In 
particular, if the central bright spot were the wind neb- 
ula (PWN) around a very young pulsar, then one might 
expect it to have a somewhat flatter spectral index, as 
mos t PWNe have spectral i ndices in the range —0.3 to 
0.0 (jGaensler fc Slanel 12006) and a s there are no filled- 
center remnants in the catalog of [Green (2004) which 
have broadband radio spectra steepeiFl than a = —0.3. 
In particular, the youngest known PWNe, the Crab Neb- 
ula and G21.5— 0.5, have r ather flat radio spectra with 
spectra l indices of —0.27 j Bictcnh olz et al.l 119971 ) and 
+0.08 (|Bietenholz fc Bartell l2008bf) . respectively. The 
observed steep spectrum of the central bright spot, there- 
fore, would not seem to favor the hypothesis of a young 
PWN. However, an only 27-yr old PWN may not be 
comparable to the much older known PWNe, so the ob- 
served steep spectrum does not rule out a young PWN. 

Alternatively, if the central bright spot represents emis- 
sion from jets emanating from a black hole environment, 
the radio spectrum is consistent with the expectations 
but the radio luminosity is not. The possible presence of 
both self-absorbed and optically thin components in the 
jet can give rise to both flat and steep spectra, so the ob- 
served spectrum can easily be accommodated. However, 
SN 1986J's central bright spot is far more radio-luminous 
than any known stellar-mass black hole system. In par- 



1 see http://www.mrao.cam.ac.uk/surveys/snrs for an up- 
dated version of the catalog. The possible exception is Vela X, 
the central source in the Vela supernova remnant, which is often 
associated with the PWN. I Alvarez etaLl pOOlT ) found a* 09GHz = 
—0.39 ± 0.03 for Vela X, indicatin g a relatively steep spectrum for 
a PWN. More recently, however, Hales et al. (2004) reported on 
31 GHz observations of a strong source near the Vela pulsar that 
they identify with the PWN, and for which they find a slightly 
inverted spectrum with «| 1 |g h z z = +0.10 ± 0.06 



ticular, a "fundamental plane" relationship between the 
radio and X-ray luminosities and black-hole mass has 
been observed for b lack- hole systems of a wide range 
of masses (see, e.g.. iFalcke ^ eTaLl[200l iHol [2008] ). The 
unabsorbed X-ray f l ux fro m SN 1986 J was measured in 
late 2003 bv lHouckl (fMOfih to be - 1.6 x 10~ 13 erg s _1 . 
Even with the assumption that all this X-ray flux is due 
to the central bright spot, the "fundamental plane" re- 
lation would suggest radio luminosities for a black-hole 
system of ~ 10~ 4 of that observed for the central bright 
spot of ~ 2 x 10 36 erg s -1 (§ 13. 3|) . However, very little is 
known about such young black-hole systems, so it is not 
impossible that they would be far more luminous than 
expected from the fundamental plane relationship. So 
again, we must conclude that the radio properties of the 
central bright spot do not seem to argue strongly for an 
accreting black-hole system, but neither can they exclude 
it. 

4.7. The Central Bright Spot as a Shell Component? 

As we showed in §[43] above, the northeast bright spot 
can naturally be explained as the impact of the shock- 
front on a dense condensation in the CSM. Could the 
central bright spot be a similar phenomenon, with the 
dense condensation being coincidentally located near the 
projected center of the shell? In this case, an approxi- 
mate equality of the intrinsic spectral indices of the cen- 
tral bright spot with that of the remainder of the super- 
nova is expected, since usually not much variation of the 
spectral index around the shell is observed. 

As we showed, there is considerable evidence that the 
CSM of red supergiants and of SN 1986J in particular are 
quite clumpy, and the northeast bright spot suggests that 
there is at least one dense condensation in the CSM of 
SN 1986J. There is no particular reason to only expect a 
single dense clump in the CSM, so the existence of more 
than one dense clump seems not unlikely. 

The expanding shell impacting on a dense CSM clump 
would cause a local brightening of the radio emission. Al- 
though the bulk of the CSM is optically thin at cm wave- 
lengths, clumps with densities as high as N e > 10 6 cm~ 3 , 
and sizes small compared to the shell diameter, would 
be expected to have large optical depths for cm- wave ra- 
dio emission. For example, a clump with a diameter of 
5xl0 16 cm, corresponding to an angle of ~0.3 mas on 
our images, and a density of lxlO 6 cm~ 3 would have a 
free-free optical depth at 8.4 GHz of ~30. 

The shock hitting a dense clump will produce localized, 
bright radio emission, as is seen in the northeast bright 
spot. However, if the dense CSM clump lies on the near 
side of the expanding shell, fortuitously near the center 
in projection, then the optically thick clump itself will 
in fact block most of this bright radio emission. With 
time, the clump will become optically thin as it either 
fragments or as the shock eats through it. This scenario 
would produce a sequence much like what is seen for the 
central bright spot: a delayed turn on, dependent on the 
distance of the clump from the explosion center and its 
density, with an inverted spectrum expected from free- 
free absorption, followed eventually by an optically thin 
decay. At 5 GHz, the brightness of the central bright spot 
is still increasing, (Fig. [7]) so we are still in the rising part 
of its lightcurve. 

If we interpret the central bright spot as emission from 
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the shock interacting with such a CSM clump, what can 
be deduced about the physical conditions in this clump? 
In our latest image at 22 GHz, we can estimate an an- 
gular diameter of the central bright spot of ^0.7 mas 
(FWHM), suggesting a linear diameter of ~ 1 x 10 17 cm 
or 0.03 pc. The turnover frequency at this epoch is 
~12 GHz. If we take the turnover frequency to be that at 
which the free-free optical depth is unity, the density of 
the absorbing clump can be estimated from the electron 
number density which in turn can be calculated from 
the emission measure using the equation given in § 14.61 
above. If we again assume T e — 10 4 K and a constant N e , 
and further assume that the clump's line-of-sight depth is 
equal to the central bright spot's diameter in the plane 
of the sky, we can calculate that N e ~ 2.5xl0 5 cm" 3 . 
Again, the presence of significantly higher densities can- 
not be excluded, but the presence of substantial mate- 
rial at the calculated density seems required. The total 
mass of such an absorbing clump would be ~0.1 M Q (as- 
suming full ionization). The swept- up mass in 2008.8 
i s 0.8 ~ 2.0 M Q (out to a radius of 6.3 x 10 17 cm; see 
§03]), so the mass of the clump represents 13% to 5% of 
the total swept-up mass. 

By comparing the part of the spectrum above 20 GHz 
with that at low frequencies (see Fig. [5]), we can esti- 
mate that the unabsorbed radio luminosity of the central 
bright spot is about twice that of the remainder of the 
supernova. Note that the unabsorbed luminosity of the 
central bright spot is currently larger than that of the 
northeast bright spot, since a substantial fraction of the 
central bright spot's luminosity is still absorbed, with 
the turnover frequency currently being ~12 GHz. If we 
again assume that radio luminosity per unit shock area 
scales only with the CSM density, we can estimate the 
density of the clump relative to that of the average CSM 
density as we did for the northeast bright spot (§ I4.5|) . 
Using the above clump diameter of lxlO 17 cm, and a 
shell outer radius of 6xl0 17 cm (for 2006 Dec.) we find 
that the spot covers ^0.2% of the shell surface. The 
clump must therefore be ~ 1000 x denser than the av- 
erage corresponding CSM. As might be expected, since 
the central bright spot has a much higher unabsorbed 
luminosity than the northeast bright spot, the density 
required for the relevant enhancement of the emissivity 
is also much higher. 

We calculated a representative density of N e ~ 
2.5 xlO 5 cm -3 for the central clump above, which sug- 
gests a representative average density for the CSM at 
r ~ 6.2 xlO 17 cm of ~ 250 cm" 3 , which in turn implies a 
mass-loss rate of 6xl0 -5 M© yr _1 (for v w = 10 km s" 1 ) 
which is consistent with the values derived in § l4.4l above. 

As mentioned in § 14.21 above, the oxygen forbidden-line 
emission seems to be dominated by two dense clumps, 
both on the near side of the remnant. If the central ra- 
dio bright spot is due to a CSM clump on the near side of 
the SN, then also in the radio, there would be evidence 
for two dense clumps. Could the clumps responsible for 
radio bright spots be be the same as the two clumps seen 
in the optical emission lines? The association seems un- 
clear, since the radio emission, as we have argued, would 
be due to dense clumps in the CSM, exterior to the for- 
ward shock, while the forbidden-line emission is plausibly 
attributed ejecta heated by the reverse shock, that is in- 



terior to the reverse shock. 

While it is conceivable that dense CSM clumps cause 
distortions first in the forward shock but subsequently 
also in the reverse shock, which latter could in turn give 
rise to the bright features seen in the oxygen forbidden 
lines, any association between the radio bright spots and 
the oxygen-line features must remain speculative. 

All said, we think now that the spectral evolution seen 
in the new observations make an interpretation of the 
central bright spot in terms of a second shell component 
and an interpretation in terms of emission from a PWN 
or a black- hole environment equally plausible. 

5. SUMMARY AND CONCLUSIONS 

1. We have obtained new multi- frequency VLA flux den- 
sity measurements and VLBI images of SN 1986J, show- 
ing the evolution of this supernova in the radio. 

2. The evolution of the integrated radio spectrum is 
complex. The spectrum at the lowest frequencies has 
a spectral index in the range of —0.7 ~ —0.5. The spec- 
trum above the high-frequency turnover, presumably the 
intrinsic spectrum of the central bright spot, is equally 
steep within the uncertainties. 

3. The shell continues to expand. The average ex- 
pansion speed of the shell between 1999 and 2009 was 
5700 ± 1000 km s _1 . This speed is compatible with con- 
tinued power-law expansion, with the radius increasing 
x increase in the rate of flux density de- 
cay at t ~ 7 yr is likely due to a flattening in the profile 
of the ejecta profile rather than a steepening in the one 
of the CSM, as no increase in deceleration is observed. 

4. The equipartition magnetic field decreases as the su- 
pernova expands, with an approximate value of B ~ 
60 (t/10 yr)" 13 mG. The dependence on the outer shock 
front radius is B cx r . 

5. Various estimates of the mass-loss converge to val- 
ues in the range of (4 ~ 10) x 10" 5 Af© yr" 1 (for v w = 
10 km s" 1 ). 

6. The VLBI images show two bright spots in addition 
to the shell-like structure. The first is in the northeast 
of the shell, while the second is near the projected center 
of the shell. The flux densities of the spots relative to 
that of the SN as a whole vary with time, with that 
of the northeast bright spot increasing till ^1999, and 
decreasing since, while that of the central bright spot 
continues to increase. 

7. The northeast bright spot is likely due to a 
dense clump in the circumstellar material (CSM). The 
bright spot's proper motion is consistent with homolo- 
gous power-law expansion. Number densities of N e ~ 
10 4 cm -3 or higher are suggested for the clump. 

8. The central bright spot has a partly absorbed radio 
spectrum, with an intrinsic radio spectrum that is sim- 
ilar to that of the shell. The amount of absorption is 
decreasing with time. 

9. The central bright spot's original interpretation as 
originating in the physical center of the shell and being 
emission due to the neutron star or black hole remnant 
of the supernova is supported by its central position, and 
its stationarity to within la. The new observations, how- 
ever, suggest an equally plausible alternative explanation 
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of the central bright spot being radio emission due to the 
shell impacting upon a second dense CSM clump, for- 
tuitously located on the near side of the shell close to 
the projected center. The amount of absorption suggests 
number densities of N e > 2.5xl0 5 cm" 3 ,~ 10 3 higher 
than the average density of the CSM. 

10. From the VLBI images from 1987 to 2008, we pro- 
duced a movie showing the supernova's evolution. 
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